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O ■ ABSTRACT 

(— I ' The degeneracy among the disk, bulge and halo contributions to galaxy rotation curves 

O . prevents an understanding of the distribution of baryons and dark matter in disk 

galaxies. In an attempt to break this degeneracy, we present an analysis of the strong 
gravitational lens SDSS J2141— 0001, discovered as part of the SLACS survey. The lens 
galaxy is a high inclination, disk dominated system. We present new Hubble Space 
^ ' Telescope multicolor imaging, gas and stellar kinematics data derived from long-slit 

spectroscopy, and K-band laser guide star adaptive optics imaging, both from the Keck 
telescopes. We model the galaxy as a sum of concentric axisymmetric bulge, disk and 
halo components and infer the contribution of each component, using information 
from gravitational lensing and gas kinematics. This analysis yields a best-fitting total 
£Tc ■ (disk plus bulge) stellar mass of log lo (M*/M0) = 10.991q'2 5 - The photometric data 

combined with stellar population synthesis models yield iog lo (M*/M ) = 10.97 ± 
0.07, and 11.21 ± 0.07 for the Chabrier and Salpeter IMFs, respectively. Assuming 
no cold gas, a Salpeter IMF is marginally disfavored, with a Bayes factor of 2.7. 
Accounting for the expected gas fraction of ~ 20% reduces the lensing plus kinematics 
stellar mass by 0.10 ± 0.05 dex, resulting in a Bayes factor of 11.9 in favor of a 
Chabrier IMF. The dark matter halo is roughly spherical, with minor to major axis 
ratio g3, h = O^I^q ^. The dark matter halo has a maximum circular velocity of 
Vmax = 276^18 km s _1 , ano - a central density parameter of log 10 A v / 2 — 5.9^q%. This 
k> , is higher than predicted for uncontracted dark matter haloes in ACDM cosmologies, 

j_j ■ logio Ay/2 = 5.2, suggesting that either the halo has contracted in response to galaxy 

formation, or that the halo has a higher than average concentration. Larger samples 
of spiral galaxy strong gravitational lenses are needed in order to distinguish between 
these two possibilities. At 2.2 disk scale lengths the dark matter fraction is /dm = 
0.55l£;^, suggesting that SDSS J2141-0001 is sub-maximal. 
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1 INTRODUCTION 

The discovery of extended flat rotation curves in the outer 
parts of disk galaxies three decades ago (Bosnia 1978; Rubin 
et al. 1978) was decisive in ushering in the paradigm shift 
that led to the now standard cosmological model dominated 
by cold dark matter (CDM). The need for dark matter on 
cosmological scales is also firmly established from observa- 
tions of the Cosmic Microwave Background, type la Super- 
novae, weak lensing, and galaxy clustering (see, e.g., Spergel 
et al. 2007). Numerical simulations of structure formation 
within the ACDM cosmology make firm predictions for the 
structure and mass function of dark matter haloes in the 
absence of baryons (e.g., Navarro, Frenk, & White 1997; 
Bullock et al. 2001; Maccio et al. 2007; Navarro et al. 2010). 

It is still unclear, however, whether this standard model 
can reproduce the observed properties of the Universe at 
galactic and sub-galactic scales. There are problems related 
to the inner density profiles of dark matter haloes (e.g., de 
Blok et al. 2001; Swaters et al. 2003; Newman et al. 2009), 
reproducing the zero point of the Tully-Fisher relation (e.g., 
Mo & Mao 2000; Dutton et al. 2007), and the amount of 
small-scale substructure (e.g., Klypin et al. 1999; Moore et 
al. 1999; Stewart et al. 2008). There are three classes of solu- 
tions to these problems: those that invoke galaxy formation 
processes that modify the properties of dark matter haloes; 
those that change the nature of dark matter itself; and those 
in which dark matter does not exist. Thus, measuring the 
density profiles of the dark matter haloes of galaxies of all 
types is a stringent test for galaxy formation theories. 

From an observational point of view, little is known 
about the detailed distribution of dark matter in the inner 
regions of disk galaxies, despite the great investment of tele- 
scope time and high quality measurements of hundreds of 
rotation curves (e.g., Carignan & Freeman 1985; Begeman 
1987; Courteau 1997; de Blok & McGaugh 1997; Verheijen 
1997; Swaters 1999; de Blok et al. 2001; Swaters et al. 2003; 
Blais-Ouellette et al. 2004; Simon et al. 2005; Noordermeer 
et al. 2005; Simon et al. 2005; Chemin et al. 2006; Kuzio de 
Naray 2006; de Blok et al. 2008; Dicaire et al. 2008; Epinat 
et al. 2008). The fundamental reason is the so-called disk- 
halo degeneracy: mass models with either maximal or min- 
imal baryonic components fit the rotation curves equally 
well, leaving the structure of the dark matter halo poorly 
constrained by the kinematic data alone (e.g., van Albada 
& Sancisi 1986; van den Bosch & Swaters 2001; Dutton 
et al. 2005). Stellar population models are able to place 
constraints on stellar mass-to- light ratios, allowing inference 
about the baryonic contribution to the overall mass profile. 
However, there are a number of uncertainties which limit 
the accuracy of this method (e.g., Conroy et al. 2009, 2010). 
These include systematic uncertainties such as the unknown 
stellar initial mass function (IMF), and the treatment of 
the various stellar evolutionary phases in stellar population 
synthesis (SPS) models. These result in about a factor of 2 
uncertainty in the stellar masses estimated from spectral en- 
ergy distribution (SED) fitting. Moreover, for a given IMF 
and SPS model, there are uncertainties in the star forma- 
tion histories, metallicities and extinction which introduce 
(la) random errors in measurements of stellar masses for 
individual galaxies at the level of 0.15 dex (e.g., Bell & de 
Jong 2001; Auger et al. 2009, Gallazzi & Bell 2009). 



Nevertheless, galaxy colours and dynamical mass esti- 
mates have been used in combination by various authors to 
place an upper limit on the stellar mass-to-light ratio nor- 
malisation, favoring IMFs more bottom light than Salpeter 
for spiral galaxies and fast rotating low-mass elliptical galax- 
ies (Bell & de Jong 2001; Cappellari et al. 2006; de Jong & 
Bell 2007; see, however, Treu et al. 2010, Auger et al. 2010, 
and van Dokkum & Conroy 2010 for massive ellipticals). 
However, as theory generally predicts more dark matter in 
the inner regions of disk galaxies than is consistent with 
standard IMFs (Dutton et al. 2007; Dutton et al. 2010b) a 
lower limit to the stellar mass would provide a more useful 
constraint for ACDM. 

Several other methods have been used to try and mea- 
sure disk galaxy stellar masses, independent of the uncer- 
tainties in the IMF. These include: 1) vertical velocity dis- 
persions of low inclination disks (Bottema 1993; Verheijen 
et al. 2007; Bershady et al. 2010), 2) bars and spiral struc- 
ture (Weiner et al. 2001; Kranz et al. 2003), and 3) strong 
gravitational lensing by inclined disks (Mailer et al. 2000; 
Winn et al. 2003). None of these methods have thus far 
yielded conclusive results. 

An approach combining strong gravitational lensing 
plus kinematics holds great promise, because it takes advan- 
tage of the different geometries of disks and haloes, which 
results in three effects that enable the disk mass to be mea- 
sured. 1) An inclined disk will present a much higher pro- 
jected surface density than a face-on disk, with resulting 
image positions and shapes that depend on the disk mass 
fraction. 2) An edge-on disk is highly elliptical in projec- 
tion, more than expected for any realistic dark matter halo, 
with resulting total mass ellipticity depending on the disk 
mass fraction 3) Strong lensing measures mass projected 
along a cylinder (within the Einstein radius), whereas stellar 
kinematics (rotation and dispersion) measure mass enclosed 
within spheres (see Figure [TJ. For spherical mass distribu- 
tions of stars and dark matter, the ratio between the pro- 
jected mass within a cylinder of radius, r, and the enclosed 
mass within a sphere of the same radius, r, is independent of 
the relative contribution of the two mass components (left 
panel Figure [2}. Therefore, in order to break the degener- 
acy one has to assume a radial profile shape for both com- 
ponents (e.g., Treu & Koopmans 2002, 2004; Koopmans & 
Treu 2003; Koopmans et al. 2006; Treu et al. 2010; Auger et 
al. 2010). Typically this involves assuming the baryonic mass 
follows the light, and then assuming a functional form for 
the dark matter halo. However, for a disk plus halo system, 
this ratio is dependent on the relative contribution of the 
two components (right panel Figure [5J). Thus if the spherical 
and cylindrical masses can be measured accurately enough, 
the disk halo degeneracy can be broken without assuming a 
specific radial profile shape for either component. Further- 
more, strong lensing plus kinematics can place constraints 
on the 3D shape of the dark matter halo (e.g., Koopmans, de 
Bruyn, & Jackson 1998; Mailer et al. 2000) which is of inter- 
est because ACDM haloes are predicted to be non-spherical 
(e.g. Allgood et al. 2006; Bett et al. 2007; Maccio et al. 2008). 

The power of the strong gravitational lensing method 
has not yet been fully realised, primarily due to the scarcity 
of known spiral galaxy gravitational lenses. Prior to the 
SLACS Survey (Bolton et al. 2006, 2008) only a handful of 
spiral galaxy lenses with suitable inclinations to enable rota- 
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Figure 1. Illustration of the different geometries probed by strong lensing and kinematics. Strong lensing measures mass with a cylinder 
(or more generally an ellipse), whereas stellar and gas kinematics measure mass within spheres (or more generally ellipsoids). 
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Figure 2. Differences between projected (cylindrical) mass and enclosed (spherical) mass for a bulge-halo system (left) and a disk-halo 
system (right). For each system two models are shown (in red and black). The models have baryonic mass profiles (short-dashed lines, 
upper panels) with the same shape but normalizations that differ by a factor of two. The dark matter profiles (long-dashed lines, upper 
panels) have been chosen so that the total circular velocity curves are close to identical (solid lines, upper panels). For the bulge-halo 
system the ratio between projected and enclosed masses (middle panels) is independent of the relative contributions of the bulge and 
halo, which differ significantly between the two models (lower panels). However, for the disk-halo system there is a significant difference 
between the projected and enclosed masses, especially at radii smaller than the effective radius. This illustrates the potential of strong 
lensing plus kinematics to break the disk-halo degeneracy. 
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tion curve measurements were known: Q2237+0305 (Huchra 
et al. 1985; Trott & Webster 2002); B1600+434 (Jackson 
et al. 1995; Jaunsen & Hjorth 1997); PMN J2004— 1349 
(Winn et al. 2003); CXOCY J220132.8-320144 (Castander 
et al. 2006). However, most of these systems are doubly- 
imaged QSOs which provide minimal constraints on the 
projected mass density. Q2237+0305 is a quadruply-imaged 
QSO, which gives more robust constraints, but since the 
Einstein radius is small compared to the size of the galaxy, 
the lensing is mostly sensitive to the bulge mass, not the 
halo (Trott et al. 2010; van de Ven et al. 2010). 

The final SLACS lens sample (Auger et al. 2009) is com- 
prised of 98 strong galaxy-galaxy lenses, among these, 16 
have been classified morphologically as type S or SO. In- 
spired by this, we have extended the Sloan Digital Sky Sur- 
vey (hereafter SDSS, York et al. 2000) spectroscopic lens 
selection technique specifically to spiral galaxy lenses. In 
the resulting SWELLS survey (Treu et al. , in prep, re- 
ferred to hereafter as Paper I) we have assembled a larger 
sample of 20 late-type galaxy-scale gravitational lenses for 
detailed mass modelling. In this paper, the second of the 
SWELLS series, we present a detailed and self-consistent 
mass model of the spiral galaxy lens SDSS J2141— 0001 
(RA=21:41:54.67, DEC=-00:01:12.2, J2000), constrained 
by both kinematic and lensing data. As we will see, this 
galaxy is disk-dominated, with a disk inclination of ~ 80°; 
this set-up approximately maximises the projected disk mass 
while allowing an accurate rotation curve to be measured. 

The original spectroscopic observations of 
SDSSJ2141-0001 were obtained on SDSS plate 989, 
with fiber 35, on MJD 52468. The latest public SDSS- 
DR7 (Abazajian et al. 2009) Petrosian magnitudes 
(uncorrected for extinction) for the lens galaxy are 
(u,g,r,i,z) = (20.61,18.62,17.47,16.92,16.48) with errors 
(0.15,0.01,0.01,0.01,0.02). The SDSS measured redshift 
for the lens galaxy is Zd ~ 0.1380 ± 0.00015, and the 
velocity dispersion is 181 ± 14kms _1 . The spectrum also 
exhibits nebular emission lines at a background redshift of 
z s = 0.7127 (Bolton et al. 2008). With these redshifts the 
scale in the lens plane is larcsec = 2.438 kpc, while in the 
source plane it is larcsec = 7.196 kpc. 

This paper is organised as follows. In Section [5] we 
present the imaging observations of SDSS J2141— 0001 from 
Keck and the Hubble Space Telescope, and then infer the 
structure of the stellar mass distribution of the galaxy in 
the presence of its dust from these data in Section [3] With 
this information in hand we then define a three-component 
mass model for the galaxy in Section U and describe how 
it is constrained by the imaging data (although we choose 
not to use the stellar mass inferred from the SED at this 
stage). In Section [5] we describe the preparation and anal- 
ysis of the strong lensing data. In Section [6] we present the 
spectroscopic observations of SDSS J2141— 0001 from Keck. 
Then in Section[7]we present fits to the lensing and kinemat- 
ics data using three combinations: lensing only; kinematics 
only; and lensing plus kinematics. This joint analysis yields 
constraints on the stellar mass of the disk and bulge; return- 
ing to the stellar masses inferred from the stellar population 
modelling of the SED we, discuss implications of our results 
for the stellar IMF in Section E3 In Section [73] & Sec- 
tion [Tj4] we discuss our results for the density and shape of 
the dark matter halo. We conclude in Section [8] 



Table 1. Summary of Imaging Observations 
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Integration Time(s) 
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F814W 
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Keck II 


NIRC2-LGS 


K' 


2700 



Throughout, we assume a flat ACDM cosmology with 
present day matter density, O m = 0.3, and Hubble param- 
eter, Ho = 70 km s _1 Mpc. All magnitudes are given in the 
AB system. Unless otherwise stated, all parameter estimates 
are the median of the marginalised posterior PDF, and their 
uncertainties are described by the absolute difference be- 
tween the median and the 84th and 16th percentiles (such 
that the error bars enclose 68% of the posterior probability). 



2 HIGH RESOLUTION IMAGING 
OBSERVATIONS 

SDSS J2141-0001 has been imaged at ~ 0.1 arcsec FWHM 
resolution from the optical (with HST) to the NIR (with 
Keck Laser Guide Star Adaptive Optics). A summary of the 
imaging observations is given in Table [T] whilst Figure 
shows the HST and Keck images of SDSS J2141-0001. In 
this section we describe the multi-filter high-resolution imag- 
ing data obtained for the SDSS J2141— 0001 system in some 
detail. 



2.1 ACS/WFPC2 imaging from HST 

Hubble Space Telescope (HST) observations of 
SDSS J2141-0001 were obtained on June 12th 2006 
with the Advanced Camera for Surveys (ACS), and on 
April 19th 2009 with the Wide Field Planetary Cam- 
era 2 (WFPC2). The ACS observation, in the F814W 
(420s) filter, was part of the SLACS snapshot programme 
GO:10587 (PLBolton); these data have a pixel scale of 
0.05 arcsec. The WFPC2 observations, in the F450W 
(4400s) and F606W (1600s) filters, were part of the cycle 
16 supplementary programme GO:11978 (PLTreu). For the 
WFPC2 observations four sub-exposures were obtained, 
and the frames drizzled to a pixel scale of 0.05 arcsec. 

The F814W image confirmed the strong lensing nature 
of this system by showing that the background object was 
multiply imaged into a three-component arc. It revealed that 
the lens was a disk dominated galaxy, with a high inclination 
and dusty disk, and that the bulge was compact and disk 
like. The F450W and F606W images reveal that the source 
is blue. 

Models for the point spread functions (PSFs) of the 
ACS and WFPC2 data where obtained by using the program 
Tiny Tim (Krist 1995). These PSF models include the effects 
of sub-pixel dithering and drizzling and have been found 
to provide adequate models for the true PSF (e.g., Bolton 
et al. 2008; Auger et al. 2009). 
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Figure 3. Optical to near-IR, high resolution imaging of SDSS J2141— 0001. Images are 12 arcsec by 6 arcsec, North up, East is left. 
The lens is a high inclination, disk dominated, and star forming spiral galaxy. The source (~ 1.3 arcsec to the east of the galaxy center) 
appears to be multiply imaged in the optical, but is a continuous arc in the near IR. 



Table 2. Summary of bulge plus disk fits together with stellar masses derived from SED fits with a Chabrier (2003) IMF. 





q 


i?50 [arcsec] 


n F450W-K' 


F606W-K' 


F814W-K' 


K' magnitude 


log 10 (M»/M ) 


Bulge 
Disk 


0.53 ±0.02 
0.31 ±0.02 


0.26 ±0.01 
2.53 ± 0.13 


1.21 ±0.11 3.80 ±0.04 
= 1.0 3.04 ±0.12 


2.40 ± 0.03 
1.75 ±0.14 


1.44 ±0.04 
1.00 ±0.12 


17.76 ±0.28 
16.25 ±0.13 


10.26 ±0.08 
10.88 ±0.07 



2.2 NIRC2 imaging from Keck 

On August 13th 2009, we imaged SDSS J2141-0001 with 
the Laser Guide Star Adaptive Optics (LGSAO) system on 
the Keck II telescope. The tip-tilt star had an R-band mag- 
nitude of 16.2 and a separation from the science target of 
60.4 arcsec. 

The images were taken in the K'-band with the near- 
infrared camera NIRC2, in wide field format (with a 40 x 40 
arcsec field of view) . The pixel scale for this configuration is 
0.04 arcsec pix -1 . A total of 45 minutes of exposure was ob- 
tained. Individual exposures were 1 minute in duration (di- 
vided into two 30-second co-adds). A dither was executed af- 
ter every set of 5 exposures to improve sky sampling. Dithers 
were based on a four point box pattern with sides 8 arcsec. 
The laser was positioned at the center of each frame, rather 
than fixed on the central galaxy. Observing conditions dur- 
ing the run were good. 

The images were processed with the CATS reduction 
procedure described by Melbourne et al. (2005). A sky frame 
and a sky flat were created from the individual science ex- 
posures after masking out all objects. Frames were then 
flat-fielded and sky-subtracted. The images were de-warped 
to correct for known camera distortion. The frames were 
aligned by centroiding on objects in the field, and finally 
co-added to produce the final image. 



A model for the PSF was derived from observations of 
a PSF star pair, where the star used for tip-tilt correction is 
the same distance from the PSF star as the lens galaxy was 
from its tip tilt star. The star pair observations were made 
immediately following the lens observations. The PSF star 
was found to have FHWM=0.10 arcsec (2.5 pixels) and a 
Strehl ratio of 18%. 

In the K'-band the extinction of both the lensed images 
and the lens galaxy light due to dust in the lens galaxy is 
almost completely absent, revealing a ring like structure, and 
confirming the disky nature of the bulge. The background 
object appears to have been lensed into a smooth arc in 
this filter. The difference between the source structure in 
the rest-frame NIR and the rest-frame UV/optical is likely 
due to extinction from the lens galaxy artificially creating 
the appearance of three distinct images. 



3 THE STELLAR MASS DISTRIBUTION 

We begin our study of the mass distribution of 
SDSS J2141-0001 by inferring the structure of the stellar 
component from the high resolution imaging data described 
in the previous section. Following the standard approach 
(e.g., MacArthur et al. 2003) our strategy is to model the 
stellar mass distribution as an exponential disk of stars plus 
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a Sersic profile bulge, with each spatial component consist- 
ing of distinct stellar populations. We first fit the surface 
brightness data to obtain estimates of the shape and profile 
of the stellar mass density, and then normalise the two pro- 
files by fitting the bulge and disk fluxes in our 4 filters (the 
spectral energy distribution, or SED) with stellar population 
synthesis (SPS) models. 

3.1 Disk/bulge surface brightness fits 

In each band, a 2-dimensional model of the lens galaxy sur- 
face brightness was fitted to the high resolution imaging 
data. The model is composed of two elliptically-symmetric 
Sersic profile components, representing the disk and the 
bulge. 

E(x,j / ) = E exph( J R/i? ) 1/n ] (1) 

where R — *Jx 2 + y 2 jq 2 . The Sersic index n is fixed at 1 
for the exponential disk, and left free for the bulge. The 
remaining parameters for each component are the centroid 
position {x c ,y c }, scale radius Ro, the axis ratio q, and the 
orientation angle (f>. The prior probability distributions were 
all independent, with uniform priors for <f), x c , y c and q, and 
"Jeffreys" (oc 1/x) priors for Ro and n, between generous 
upper and lower bounds. 

All four bands are fitted simultaneously, with all pa- 
rameters except for the normalization of the bulge and disk 
fluxes constrained to be the same in all bands. This approach 
gives more robust colors of the bulge and disk than is ob- 
tained when letting the structural parameters float between 
bands. 

The inferred parameter values for the disk and bulge 
surface brightness are given in Table [2] For the bulge com- 
ponent, we find a Sersic index of rib = 1.21±0.11, and bulge 
(luminosity) fraction which increases from 0.11 ±0.03 in the 
F450W filter to 0.20 ±0.05 in the K' band. These values are 
typical for low-redshift late-type spiral galaxies. 

The bulge has a major-axis half-light radius of J?5o,b = 
0.26 ± 0.01" = 0.63 ± 0.02 kpc, whilst the disk has a major- 
axis half-light radius of R 30 ,d = 2.53±0.13" = 6.17±0.32kpc 
corresponding to a disk scale length Rd = 1.51" ± 0.08" = 
3.68±0.19 kpc. The ratio between the bulge half-light radius 
and the disk scale length is 0.17 ± 0.02 which is consistent 
with those found by MacArthur et al. (2003) in a sample of 
moderately inclined late- type spirals. 

The bulge has an observed axis ratio of qt = 0.53±0.02, 
while the disk has an observed axis ratio of q^ = 0.31 ±0.02. 
For a thin disk the axis ratio equals the cosine of the 
inclination angle: g<j = cos(i). However, in general disks 
have a finite thickness, which causes the true inclination to 
be higher than that inferred from the observed axis ratio. 
For SDSS J2141— 0001 we can infer the disk inclination by 
measuring the axis ratio of the star forming ring, and as- 
suming it is intrinsically circular. This yields an axis ratio 
q = 0.20 ± 0.02, and thus i = 78.5 ± 1.2 degrees. 

For an oblate ellipsoid with projected minor-to-major 
axis ratio, q, and inclination, i, the 3D minor-to-major axis 
ratio, <?3 is given by 

qf = (q 2 — cos 2 «)/(l — cos 2 i) (2) 

Thus for SDSS J2141-0001 we infer that the 3D axis ra- 
tio of the bulge and disk are q^b = 0.51 ± 0.02, and 
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Figure 4. Posterior distributions for the inference on the stellar 
mass based upon SPS models constrained by the high-resolution 
photometry, assuming a Chabrier IMF. 

q-i, d = 0.26 ± 0.02. The disk thickness that we derive for 
SDSS J2141— 0001 is in good agreement with measurements 
of edge-on spiral galaxies (Kregel et al. 2002). We note that 
the orientation angles of the disk and bulge components are 
very close to each other: SDSS J2141— 0001 appears to be 
well-modelled by an oblate bulge bisected by a thick, coax- 
ially symmetric disk. 

3.2 Stellar population SED fits 

The two component model for the surface brightness of 
SDSS J2141— 0001 inferred in the previous section can be 
used to constrain the stellar mass distribution of the galaxy. 
Given the precision of the shape and profile measurements, 
the normalisation of the stellar mass distribution is the most 
uncertain part. Our aim is to constrain this both gravita- 
tionally (via dynamical and lensing measurements) , and by 
modelling the stellar populations of the disk and bulge. In 
this section we describe the latter route. 

Assuming each spatial component has a distinct stel- 
lar population, we can fit the model photometry for each 
component to SPS models and infer the stellar mass of each 
component using the code describe in Auger et al. (2009). 
We consider stellar populations characterised by either a 
Chabrier (2003) or Salpeter (1955) IMF and described by 5 
parameters: the total stellar mass M*, the population age A, 
the exponential star formation burst timescale r, the metal- 
licity Z and the reddening due to dust, tv- We employ a 
uniform prior requiring 9 ^ logio(M*/MQ) ^ 13, the age 
is constrained such that star formation began at some (uni- 
formly likely) time between 1 ^ z ^ 5, r has an exponential 
prior with characteristic scale 1 Gyr, and we impose uniform 
priors on the logarithms of the metallicity and dust extinc- 
tion such that —4 ^ logioZ ^ —1.3 and —2 ^ logioTy ^ 0.3. 
We note that the priors are the same for the bulge and the 
disk components but are sufficiently conservative that they 
do not bias our results. The posterior PDF is sampled as 
described in Auger et al. (2009). 

Figure [4] shows the marginalised posterior inference on 
the stellar mass, which we find to be well-constrained for 
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the bulge and disk components. Assuming that both com- 
ponents are well described by a Chabrier IMF, we find 
log 10 (Af*, b /M Q ) = 10.26 ± 0.08 and log 10 (M,, d /M©) = 
10.88 ± 0.07, for the bulge and disk respectively, justifying 
our description of SDSS J2141— 0001 as "disk-dominated." 
The total stellar mass of SDSS J2141-0001 from SED fit- 
ting is therefore log 10 (M„/M©) = 10.97 ± 0.07, and the 
bulge fraction is /bulge = 0.20 ± 0.04 (the same as in the 
K'-band light). For a Salpeter IMF the masses are all 0.24 
dex higher. We will return to this inference in Section 17.21 
below, where we compare it to the stellar mass implied by 
the gravitational analysis. 



4 A THREE-COMPONENT GALAXY MASS 
MODEL 

As indicated in the previous section, it makes sense to con- 
sider a stellar mass distribution for SDSS J2141— 0001 whose 
profile and spatial shape is tightly constrained by the struc- 
tural analysis of the surface brightness observed. However, 
we would like to remain agnostic about the stellar IMF, 
and instead investigate the stellar mass of the galaxy inde- 
pendently, using strong lensing and stellar kinematics. Since 
these are sensitive to the total mass of the galaxy, we include 
a dark matter halo component to the model as well. In this 
section we describe this 3-component mass model in some 
detail, including the predictions it makes for the observable 
effects, and the prior PDFs we assign on the model param- 
eters. 

4.1 Description 

Based on the results of the previous section, we model the 
stellar mass distribution of SDSS J2141— 0001 as a thin, cir- 
cular exponential disk of stars of mass M„.d, co-axial with 
an oblate bulge of stars of mass M* t b- We then assume the 
galaxy to reside in a dark matter halo that is also axisym- 
metric, and aligned and concentric with the disk and bulge. 
This assumption that the galaxy and inner dark matter 
haloes are aligned is supported by cosmological simulations 
of disk galaxy formation (e.g., Deason et al. 2011). We note 
that for our strong lensing analysis it is feasible to allow the 
position angles of the baryons and dark matter to be offset. 
However, this would make the model non-axisymmetric and 
thus make the kinematics considerably harder to model. The 
surface brightness in our four filters constrains the spatial 
distribution of stellar mass tightly, under the assumption 
that the stellar mass-to-light ratios are radially constant; 
we leave the overall normalisation of the stellar mass distri- 
bution as a free parameter. 

We do not explicitly include a cold gas disk for two 
reasons. Firstly, we do not have direct observations of the 
atomic and molecular gas in SDSS J2141— 0001. Secondly, 
observations suggest that the gas fractions for low-redshift 
luminous spiral galaxies are of order 20% (e.g., Dutton & van 
den Bosch 2009), and thus do not contribute significantly 
to the baryonic mass. However, if we were to assume that 
any cold gas present has the same spatial distribution as 
the stars, the unknown gas mass could be absorbed into a 
baryonic mass-to-light ratio that includes the stellar mass- 
to-light ratio. For the majority of this paper we neglect the 



cold gas mass, but return to it in the discussion of the IMF 
(Section [73]) below. 

The three mass components are described as follows: 

• Exponential Stellar Disk 

• Sersic Stellar Bulge 

• Non-singular Isothermal Ellipsoid (NIE) Dark Matter 
Halo 

This model has 17 parameters in total; they, and their prior 
PDFs, are given in Table [3] We assign informative priors to 
all but 4 of these parameters, propagating the uncertainties 
in the surface brightness fits through to the mass model. 

First, we assume the bulge, disk, and halo inclination 
are all the same, and given by the thin disk axis ratio (0.2), 
as in Section [3] - we assume that this is known with no un- 
certainty. As we describe below, we use an approximation to 
the Exponential profile that allows us to compute predicted 
observable quantities efficiently - the size parameters of the 
bulge and disk in that approximation are determined from 
the results of the previous section, as is (more straightfor- 
wardly) the bulge axis ratio. We use the derived value of 1.21 
(see Table[2j for the Bulge Sersic index, with no uncertainty. 

We assume that the disk and bulge are different stellar 
populations, and so use the independent stellar mass results 
from the previous section to constrain the bulge mass frac- 
tion, /bulge = M t ^/M*. As already mentioned, we leave the 
total stellar mass M, as a free parameter with uniform prior 
on its logarithm. This is effectively equivalent to assuming 
that the two components have very similar, although un- 
known, IMF normalization, We do inform the bounds of this 
uniform prior using the SPS modelling results, in the follow- 
ing way. Estimating that the lightest conceivable IMF would 
give stellar masses systematically a factor of two lower than 
Chabrier, we take the 3-sigma point of the Chabrier PDF 
in Figure [4] and subtract 0.3 dex to set a lower limit on 
log 10 (M./Mq) of 10.5. Likewise, at the high end we take 
Salpeter to be the heaviest IMF and use the 3-sigma point 
of the Salpeter PDF in Figure [4] (11.4) as our upper limit 
on log 10 (M,/Mq). We note that none of our results change 
if we adopt a higher upper limit to the stellar mass. 

This leaves 3 parameters that describe the model dark 
matter halo: V c ,h, r Ct h, and q3,h- We allow the axis ratio of 
the halo to be greater than unity, corresponding to a pro- 
late halo, but use a broad lognormal distribution centred on 
spherical to encode approximately our expectations. For the 
halo density profile, the NIE profile has considerable free- 
dom and can represent a much broader range of behaviours 
than those seen in simulation. Therefore, we adopt physi- 
cally motivated priors to select the cosmologically motivated 
subset of parameters combination. Studies of large sets of 
spiral galaxies, using satellite kinematics and weak galaxy- 
galaxy lensing, in the context of numerical simulations have 
shown that the maximum observed circular velocity is typ- 
ically comparable to the maximum circular velocity of the 
halo, even though these two maximums occur at vastly dif- 
ferent radii (Dutton et al. 2010a). We also know that ro- 
tation curves do not keep rising indefinitely, but typically 
flatten out within a few scale radii of the disk. To inject this 
information we require that the asymptotic circular velocity 
of the halo be comparable to that measured via spectroscopy 
(see §[6] below): in practice we assign a broad Gaussian prior 
centred on 280 kms -1 with width 50 kms -1 . The prior is 
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Table 3. Summary of mass model priors. Bold indicates an uninformative prior, 
regular text indicates an informative prior (i.e. that the parameter is virtually 
fixed). LAf(a,b 2 ) denotes a lognormal distribution, with a being the central value 
for the variable, and b being the standard deviation for the log of the variable. 
Af(a, b 2 ) denotes a normal distribution, with a being the central value and b being 
the standard deviation. U (a, b) denotes a uniform distribution with lower and 
upper limits, a and b, respectively. For clarity we have arranged the parameters 
into three groups: free parameters first (with uninformative priors), stellar mass 
distribution parameters (with priors derived from the SPS modelling of Section |3j, 
and finally nuisance parameters. 



parameter 


description 


prior 


Vr >, / km s — 1 


dark halo asymptotic circular velocity 


Af(280, 50 2 ) 


93, h 


dark halo 3D axis ratio 


UV(1, 0.3 2 ) 


r c h/arcsec 


dark halo core radius 


U(0.01, 10) 


log 10 (M*/M ) 


stellar mass 


U(10.5,11.4) 


/bulge 


bulge stellar mass fraction 


Af(0.2,0.04 2 ) 


9bulgc 


bulge 2D axis ratio 


ZJV(0.53,0.03 2 ) 


-Ro,bul gc /arcsec 


bulge chameleon size 


ZJV(0.094,0.03 2 ) 


"bulge 


bulge chameleon index 


0.4892 


<7disk 


disk 2D axis ratio 


UV(0.31,0.03 2 ) 


^0,disk/arcsec 


disk chameleon size 


ZjV(1.10,0.03 2 ) 


"disk 


disk chameleon index 


0.63 


cos(i) 


cosine of disk inclination angle 


0.2 


x c /arcsec 


spatial offset in x direction 


Af(0,0.01 2 ) 


j/ c /arcsec 


spatial offset in y direction 


A/"(0,0.03 2 ) 


6»/deg 


mass-light position angle offset 


Af(1.7, 2.9 2 ) 


7 


lens external shear 


A/"(0,0.1 2 ) 



# 7 /deg position angle of external shear t/(0,180) 



In later sections we will introduce the kinematic and 
lensing data, and then use them to constrain the parameters 
of this 3-component mass model. However, before getting to 
the data, in the rest of this section we give the functional 
forms for each mass component, and the predicted observ- 
ables resulting from them. 



chosen to be broad enough - the 3-sigma range of this Gaus- 
sian spans the range 130 to 430 kms" - not to drive the 
final inference and yet tight enough to rule out models where 
the maximum velocity is reached too far out. In addition we 
impose a uniform prior PDF for the core radius, allowing it 




angles can be computed analytically (Keeton & Kochanek 
1998). This model has been used very successfully to model 
the total (dark plus stellar) mass profiles of elliptical galaxy 
lenses (e.g., Bolton et al. 2008). In this work we use the NIE 
model for the halo alone. While the NIE profile has a con- 
stant central density, it is flexible enough to broadly capture 
the change in the density profile in the central regions that 
we expect from numerical simulations of dark matter halos 
(e.g., Navarro et al. 1997). 



The axisymmetric ellipsoidal halo is assumed to have a non- 
singular isothermal (NIE) profile, which we parametrise in 
a cylindrical coordinate system in the plane of the galaxy 
following Keeton & Kochanek (1998): 



4.2 Three-dimensional component density profiles 




We would like to model the stellar disk and stellar 
bulge mass components such that in projection they appear 
to have exponential and Sersic profiles respectively. How- 
ever, we also need 3D distributions for which we can com- 
pute predicted rotation curves, as well as projected distri- 
butions convenient for lensing calculations. To achieve this 
we note that the NIE profile can be used to create an ap- 
proximation to an exponential profile in projection (Mailer 
et al. 2000). This is done by taking the difference of two 
NIEs. If pnie(-R, z; V c ,r c ,q3) is a softened isothermal ellip- 
soid, then 



Here, V c is the asymptotic circular velocity, r c is the core 
radius, 53 is the three dimensional axis ratio, and e = 
(1 — q'i) 1 /' 2 is the eccentricity. For a zero thickness mass 
distribution (53 = 0), e/ sin -1 e = 2/n. For a spherical mass 
distribution (93 = 1), e/ sin -1 e = 1. For a prolate mass 
distribution (§3 > 1), while e is imaginary, e/ sin -1 e is real 
and greater than 1. 





Pnie(-R, z; V c ,r c /a,q3) 



This mass profile is often used in gravitational lens 
analysis, since its projected mass distribution and deflection 



is a "Chameleon" profile with positive density everywhere, 
and a finite total mass. In Appendix[A]we derive new formu- 
lae that provide Chameleon approximations to Sersic profiles 
of any index (for 1 < n < 4), providing the Chameleon size 
r c and index a given a Sersic half light radius R50 and in- 
dex n. 
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4.3 Predicted rotation curves 

For our ellipsoidal mass profiles, we can calculate the rota- 
tion velocity, as a function of radius, of a massless test parti- 
cle moving on a circular orbit in the plane of the galaxy. We 
refer to this velocity as the circular velocity to distinguish it 
from the rotation velocity of the stars and gas, which may be 
lower than the circular velocity due to a velocity dispersion 
component. The circular velocity profile for the NIE model 
is (Keeton & Kochanek 1998) 



V^ 1E (R;V c ,r Cl q 3 ) 
V 2 



x tan 



sin^e (R 2 + e 2 r 2 y/ 2 
(R 2 + e 2 r 2 )^ 2 



(5) 



where again e = (1 — q 2 ) 1 / 2 is the eccentricity of the mass 
distribution and the model is normalised so that, asymptot- 
ically for R — > oo, Vnie{R) — > V c . For the special case of a 
zero thickness mass distribution (q 3 = 0, e = 1) Equation [5] 
reduces to 



V&E(g) 

v 2 



(R 2 +r 2 ) 1 / 2 ' 



(6) 



For the case of a prolate mass distribution (g3 > 1), 
e/sin _1 e is real, but since e 2 < 0, for (R 2 + e 2 r 2 ) < 
the circular velocity is given by 



V 2 



= 1 



sinh 



x tanh 



e (-R 2 
(-R 2 - , 



e 2 r 2)l/2 



<73 r c 



(7) 



where e = yje^J. 

For the chameleon profile the circular velocity is given 
by the quadratic difference between the circular velocities of 
the sub-component NIE's: 



V chsn (R;V c ,r c ,q 3 ,a) = V mE (R;V c ,r c ,q 3 ) 
-V^ IE (R;V c ,r c /a,q 3 ). 



(8) 



Likewise, for the mass model the total circular velocity is 
given by the quadratic sum of the circular velocities of the 
bulge, disk, and halo components: 

V 2 {R) = V^ se (R) + v£ Bk (R) + V^ lo (R) (9) 



4.4 Predicted lensed images 

Projecting the three components onto the sky allows us to 
compute deflection angles and predict the observed gravita- 
tional arc, pixel by pixel. 

In projection the mass distribution (an oblate or prolate 
ellipsoid with minor to major axis ratio q 3 ) has projected 
axis ratio q given by 



(q 3 sin 1 + cos i, 



1/2 



(10) 



where i is the inclination angle (such that i = 0° corresponds 
to a face-on disk, and i = 90° to an edge-on one). In general 
the projected axis ratio, q, will be closer to unity than the 
3D axis ratio, q 3 . 



The projected mass density of an NIE model is given 
by (Keeton & Kochanek 1998): 



£nie(3;,2/;&, r c ,q) 



V 2 



AGD d sin" 1 e 

£crit& 



q^Jr 2 + x 2 + (y/q) 2 ' 

(11) 



2 q\Jr 2 + x 2 + (y/q) 2 

Here, Dd is the angular diameter distance to the lens, and e 
is again the ellipticity, while in the second line 6 is the minor 
axis of the critical curve (and thus b/q is the major axis of 
the critical curve), and £ cr j t is the critical surface density of 
strong lensing: 

2 



(12) 



4ttG D is D d ' 

where D s is the angular diameter distance from the observer 
to the source, and Dd s is the angular diameter distance 
from the lens to the source. For our assumed cosmology, 
for SDSS J2141-0001 these distances are: D d = 497.6Mpc, 
Ds = 1510. 2Mpc, Dds = 1179. 6Mpc, and thus the critical 
density is E crit = 4285.3M pc -2 . 

To explain the parts of Equation [TT] a little further, the 
parameter, b, is related to the spherical Einstein radius, 6sis, 
via: 



6sis(e/ sin 



(13) 



and the spherical Einstein radius (in radians) is in turn re- 
lated to the asymptotic circular velocity, V c , via: 



ftsis = 2n(V c /cYD ds /D s 



(14) 



The deflection angles are given by (Keeton & Kochanek 



1998) 



(l-g2)l/2 

b 

(1 - q 2 ) 1 / 2 



tan 



tanh 



(1 



* + r- c 

(i- q 2 ) 1/2 y 



* + q 2 r c 



(15) 

(16) 
(17) 



where ty 2 = q 2 (r 2 + x 2 ) + y 2 . The deflection angles from 
the three components of the mass distribution can be sim- 
ply summed, as they are just the first derivatives of the 
projected (lens) potential of each component, and the po- 
tentials of the three components can be summed themselves 
to give the total potential. Likewise, the Chameleon profile 
in projection is just the difference between two projected 
NIE models, and its deflection angles are just the difference 
between those of its NIE components. 

To predict the positions and structure of the lensed im- 
ages given a set of mass model parameters, we map each 
observed pixel location back to the source plane using the 
overall deflection angle map, and look up the surface bright- 
ness of a model source at that position. In practice we use a 
single, elliptically symmetric source with a Sersic brightness 
profile, as in e.g., Marshall et al. (2007). 



5 STRONG GRAVITATIONAL LENSING 
DATA 

We now present the strong gravitational lensing data that 
we will use to constrain our mass model. We first describe 
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Modified Lensing Data 



-< >• 

5" 

Figure 5. Upper panel: Galaxy subtracted image in the K'-band 
obtained using a reflection of the image. The lens galaxy is well 
subtracted near the arc, but there are significant residuals near 
the center. Lower panel: An arc in the right place, obtained by 
masking out non-arc features from the image in the upper panel. 
This was used for the actual fitting, to weaken the fit criterion. 
Essentially, we want to produce models whose posterior distribu- 
tion predicts an arc with this morphology, but does not necessarily 
need to match every pixel. 

the preparation of the arc imaging data, and then show with 
a simple lens model the information it contains. 

5.1 The lensed arc likelihood 

Due to the strong effects of dust in the lens system, we focus 
our lensing analysis primarily on the K'-band NIRC2 image. 
In the K'-band the lens galaxy appears to be much smoother 
than in the optical, but the light distribution is not able 



to be modelled by a simple surface brightness profile. This 
makes subtracting the lens galaxy light difficult. Our goal 
is to obtain robust parameter inferences with meaningful 
uncertainties, and so we opt for quite a conservative appli- 
cation of the imaging data. To account for lens subtraction 
errors, we create an arc image and a goodness of fit statistic 
that rewards a model for having an arc in the right place, 
with the right shape, and that is all: we do not require the 
detailed features of the modified surface brightness profile 
to be matched. 

To achieve this, we first subtracted the galaxy light 
around the arc by reflecting the galaxy along the minor 
axis. This method provides a better subtraction than mul- 
tiple Sersic components, or a radial bspline model (as used 
by e.g., Bolton et al. 2006). We then cut out the arc and set 
the remaining pixels to zero, before adding noise at the level 
of <ri5 = 15% of the peak arc brightness. This 15% value is 
an initial estimate of the appropriate noise level needed to 
suppress models that predict significant lensed features else- 
where, although faint counter-images are still allowed. The 
resulting modified image is shown in Figure [5] (lower panel); 
we also show the lens subtracted image (upper panel), with 
its uncertain central region. It is not clear whether or not 
there is a counter-image in the centre of the system. We note 
that if there is no counter image, then SDSS J2141— 0001 
would be a rare example of galaxy-scale naked cusp lens. 

The likelihood function for the modified image data was 
that used by Brewer & Lewis (2006), Marshall et al. (2007), 
and others. We assume Gaussian errors of 015 on the pixel 
values d; we can predict these pixel values from a model 
source with parameters S (as described in Section 14.41 
above) given lens model parameters # m . Denoting the pre- 
dicted data as d p , we write down the usual chi-squared misfit 
function 

* - 1^ 2^ ■ (18) 

We allow the data to inform our understanding of the model 
uncertainty, by re-scaling the denominator by a factor T. 
This corresponds to increasing or decreasing the perceived 
errors on the pixel values, and provides a mechanism for 
avoiding over-fitting the arc structure or allowing models 
that predict undetected flux. (The symbol T stands for 
"temperature" - increasing the temperature increases the 
diffusion of the model around its parameter space.) The like- 
lihood function is then: 

Pr(d|0 m ,0 s )(xexp(-^ X 2 ) (19) 

The value of T selected was 7.5. This was the highest 
value where the posterior distribution only contained im- 
ages that resemble the arc morphology. Higher temperatures 
caused the posterior models to predict substantial flux that 
is not observed, lower temperatures enforced the fit to the 
modified image to be too strict. The reason for the two-step 
procedure (adding noise and then selecting a temperature) 
is that Nested Sampling provides the results for all temper- 
atures in a single run, whereas tuning the noise level itself 
would have required large numbers of trial runs. 

Using this modified image, and the temperature-raising 
scheme, allows us to explore an approximate posterior dis- 
tribution for the lens model parameters that conditions on 
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Figure 6. Marginalised posterior PDF for the 2D axis ratio, q, 
and circular velocity, V c , for a single SIE component model fitted 
to the strong lensing data. The dashed lines show the median 
values of q and V c . 

the presence of an arc with the observed morphology, and 
nothing else. 

5.2 Constraints on projected ellipticity and mass 

To illustrate the unique information that strong gravita- 
tional lensing provides, we first perform a fit to the lensing 
data with a single singular isothermal ellipsoid (SIE) mass 
model. The purpose of this exercise is to show that strong 
lensing places constraints on the axis ratio of the projected 
mass, as well as the projected mass within the Einstein ra- 
dius. 

Fixing its centroid and orientation to that of the lens 
surface brightness, our example SIE lens model has two pa- 
rameters, minor axis Einstein radius 6, and axis ratio q. 
We assign uniform priors over wide ranges for these pa- 
rameters, and then explore the posterior PDF using our 
sampling code (which we introduce in more detail in Sec- 
tion [7] below) . We find the circularised Einstein radius to 
be #Ein = b/y/q — 0.891q'q| arcsec, and the axis ratio to 
be q = 0.42^0 i2- We can transform samples in b and q 
into the circular velocity, which we expect to be well con- 
strained. Figure [6] shows the marginalised posterior PDF for 
q and 14 - the circular velocity is indeed well-constrained: 
V c = 254^8 kms" . The shape of the arc also constrains the 
ellipticity of the total mass distribution: since in our three- 
component model the ellipticity of the disk and bulge are 
fixed, we expect strong lensing to then provide information 
about the shape of the dark halo. 

6 GAS AND STELLAR KINEMATICS DATA 

The second set of data that we will use to constrain our 
mass model is a galaxy rotation curve, derived from optical 



emission and absorption line spectroscopy. In this section we 
describe the observations, and the rotation curve extraction 
process, discuss the observed velocity dispersion and our in- 
terpretation of it, and then derive the likelihood function 
that we will use when fitting our mass model. 

6.1 Spectroscopic observations with Keck 

Major axis (PA = 87°) long-slit spectra were obtained with 
the DEep Imaging Multi-Object Spectrograph (DEIMOS), 
and Low Resolution Imaging Spectrograph (LRIS) on the 
Keck 10-m telescopes. 

On October 1st 2008 SDSS J2141-0001 was observed 
with DEIMOS on Keck II. We used the 1200 line grating 
(corresponding to a pixel scale of 0.32A) with a 1" width 
slit resulting in a spectral resolution of ~ 1.9 A. The central 
wavelength was 6500A, resulting in a wavelength range of 
5200 — 7800A We took three exposures of 1200s in seeing 
conditions of 0.60". The slit was aligned with the major axis 
of the galaxy, with PA = 87. The spectra were reduced using 
routines developed by D. Kelson (Kelson 2003). 

On November 27th 2008 we observed SDSS J2141-0001 
with LRIS on Keck I. However, the seeing for this ob- 
servation was considerably worse (~ 1.5") than that of 
our DEIMOS observation. This resulted in increased beam- 
smearing and reduced sensitivity, and thus we focus our 
kinematic analysis on the DEIMOS observations. 

6.2 The observed rotation curve 

Cutouts of the DEIMOS long-slit spectrum centered around 
prominent emission and absorption lines are shown in Fig- 
ure [7] These show clear signs of rotation in both emission 
and absorption lines. For the emission lines we measured 
rotation curves by locally fitting Gaussian line profiles to 
one-dimensional spectra extracted along the slit. For the ab- 
sorption lines we measured the rotation and dispersion pro- 
file by applying python routines developed by MW. Auger 
to one-dimensional spectra extracted along the slit. The up- 
per right panel also shows the spatially offset [O n] emission 
lines from the source galaxy. 

The extracted rotation curve is shown in the upper pan- 
els of Figured The spatial sampling is ~ 0.59" (5 DEIMOS 
pixels), corresponding to 1 data point per seeing FWHM. 
There is good agreement in the rotation curves measured in 
Ha and Nil, except near the very center, where Nil gives a 
higher Kot- This is possibly due to Ha being contaminated 
with stellar absorption. The Na D and Mg b absorption lines 
give lower Kot than the emission lines, especially at larger 
radii. This is expected due to the increased pressure support 
in the stars compared to the gas, so called asymmetric drift. 

The rotation curve flattens out beyond 3 arcsec (7 kpc), 
corresponding to 2 disk scale lengths. On the East side 
the rotation curve remains flat out to the last data point 
(13 kpc). On the West side the rotation curve decreases be- 
yond 3.5 arcsec (8.5 kpc). We trace this asymmetry to the 
warp in the West side of the optical disk, which causes the 
slit to miss the major axis: beyond 3.5 arcsec we therefore 
use only the East side of the rotation curve. 

Figure [9] shows the folded rotation curve obtained by 
combining the rotation curves from Ha and Nil. The data 
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Figure 7. Cutout images of the DEIMOS optical long slit spectrum of SDSS J2141— 0001 centered on prominent emission and absorption 
features: Mgb 5177, Fell 5270 from the lens (upper left); NaD 5896 from the lens (lower left); Oil 3727 from the source (upper right); 
and Ha 6563, [Nil] 6550 6585 from the lens (lower right). The vertical scale is 20 arcsec; the vertical bar is 3 arcsec. Up corresponds to 
East (left in the images in Figure [3jl- The ring in the imaging (Fig.|3jl corresponds to strong Ha emission due to star formation. 



points shown in this figure are given in Tabled When com- 
bining data points we use the error weighted mean. The 
new error is the maximum of the statistical error and half 
the differences between the two data points. From this ro- 
tation curve the maximum observed rotation velocity is 
271 ± 4kms -1 (corrected for inclination, but not beam 
smearing) at 13 kpc from the galaxy center. 

6.3 The rotation curve likelihood 

For a given set of mass model parameters G m we can predict 
the circular velocity of the stars at each radius, and hence 
find parameter vectors that fit the rotation curve data. To 
do this we need to fold in the effects of beam-smearing; our 
procedure for this is described in the next section. We com- 
pare the predicted data v p and the observed data in the 
usual way. We assume uncorrelated Gaussian errors on the 
observed velocities Vj Oj — \J K<j§,j + er^ xtra (where <To,j are 
the reported error bars, and K and <7 ex tra are free param- 
eters allowing the true uncertainties to be inferred), and 
hence construct the familiar likelihood function 



Pr(v|0 m ) = 



where the misfit function 



(20) 



(21) 



We then take the product of this likelihood and the prior 
PDFs on the parameters defined in Section [4] to obtain the 
posterior PDF for the model parameters given the kinemat- 
ics data. 
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6.4 Modelling beam-smearing 

The rotation curve data presented in Figure [9] are the ob- 
served values, uncorrected for inclination, finite slit width 
and seeing effects. We refer to these combined effects as 
beam- smearing. Since the disk inclination is high, the sini 
correction is small, just a factor of 1.021. However, since the 



Figure 8. Rotation curve (upper panel), velocity dispersion pro- 
file (middle panel), and line flux profile for SDSS J2141-0001. 
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Figure 9. Folded rotation curve from emission lines of Ha and 
[Nil]. The location of the bulge half-light radius, i?b, Einstein 
radius, -REm, and 2.2 disk scale lengths, are indicated with 
dotted lines. The spatial sampling is one point per seeing FWHM, 
which is indicated by the yellow circle. This rotation curve is 
uncorrected for inclination, and beam smearing due to the finite 
(1 arcsec) slit width and seeing. 



slit width covers a large fraction of the minor axis of the 
galaxy, the effects of finite slit width and seeing are likely 
to be significant, especially near the centre of the galaxy. 
We take this into account when computing the predicted 
data v p as the inclined, beam-smeared, model rotation curve 
within a 1 arcsec slit. For computational efficiency we esti- 
mate the beam-smearing effect using a simplified, rotating 
exponential disk model, and then apply this correction to 
the model rotation curve. The intrinsic rotation curve v is 
given by the sum of the bulge, disk and halo components as 
described above. 

The beam-smearing calculation is approximate, because 
we don't know the exact distribution of the Hcv emission, 
only the starlight. While we model the Ha distribution with 
an exponential profile, with the scale length of the V-band 
light, the actual distribution is likely to be asymmetric (due 
to extinction), and non-exponential (there is a ring of star 
formation). To minimize the impact of these uncertainties, 
we have excluded the inner 2 arcsec of data in our mass 
models. 



6.5 The central velocity dispersion 

The spectral line fits described in the previous section also 
yield some information on the velocity dispersion of the sys- 
tem. The central (within 1") velocity dispersion from the 
Mgb-[Fen] lines was found to be a = 180 ± 4kms~ 1 , 
in agreement with the SDSS value (which is integrated 
over the 3" fibre aperture). NaD gives a lower central ve- 
locity dispersion, of a = 119 ± 6kms~ 1 . Absorption in 
NaD can come from interstellar gas, as well as stars. Since 
SDSS J2141-0001 has a dusty gas disk, it is thus likely that 
the NaD line is not reliably tracing the stellar velocity dis- 
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persion. For the emission lines the central (within the inner 
0.25") velocity dispersion of the Nil line is ~ 178±14kms _1 , 
similar to that of the stars. However, the velocity disper- 
sion of the N n line declines faster with radius than that of 
the stars (middle panel of Figure [8]). This is an indication 
that the peak in velocity dispersion in Nil is due to beam- 
smearing. For the Ha line the central velocity dispersion is 
considerably lower than that of N II, which we ascribe to the 
presence of absorption, which we have not corrected for. In 
the outer part of the disk, the observed velocity dispersion 
of the emission lines is close to that of the instrumental reso- 
lution of ~ 32 kms -1 (dotted horizontal line in middle panel 
of Figure [HJ, indicating that the intrinsic velocity dispersion 
of the line-emitting gas disk is too low to be resolved. 

How could we model the velocity dispersion data? Our 
simple dynamical model does not easily allow for this, but 
we can make use of the dispersion information as a cross- 
check in the following simplistic way. The results of Pad- 
manabhan et al. (2004) and Wolf et al. (2010) show that, 
for a spherical system, the circular velocity at the half- 
light radius V c irc(^5o) — 1.7<ti os , where <ti os is the inte- 
grated line of sight velocity dispersion of the system. For 
the case of SDSS J2141-0001, the bulge half light radius is 
Rso,b — 0.26", and thus the integrated velocity dispersion 
within the inner 1" gives the integrated velocity dispersion 
of the bulge. 

Applying the Vdrc formula above and adopting an un- 
certainty of 10% results in an estimate of the circular veloc- 
ity at i?5o of Vri rc (-R5o,b) = 306 ± Slkms" 1 . In our current 
analysis we do not make use of this constraint. Rather, we 
use this as a consistency check to our models which are con- 
strained by strong lensing and gas rotation curve. 



7 RESULTS 

We now infer the parameters of our 3-component mass 
model using constraints from the strong lensing and kine- 
matics data presented in the previous two sections. In or- 
der to sample the posterior distribution for the parameters 
we use the Diffusive Nested Sampling code from Brewer 
et al. (2009). Diffusive Nested Sampling is a powerful and 
efficient alternative to standard MCMC sampling. 
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Figure 10. Marginalised two-dimensional posterior PDFs for unconstrained mass model parameters (see Table[3]l using constraints from 
strong lcnsing alone. In the histogram panels, the vertical red lines show the median value. The median value together with the offsets 
to the 84th and 16th percentiles of the distribution is given in the top right corner. The priors are shown with solid red lines. 



7.1 Inferred Model Parameters 

We consider inferences from three data sets: 

(i) strong lensing only 

(ii) kinematics only 

(iii) strong lensing plus kinematics 

In Figures [1014121 we plot, for each data set, all possi- 
ble one-dimensional and two-dimensional marginalised pos- 
terior PDFs for the four main mass model parameters. These 
parameters are the total (disk+bulge) stellar mass M* , and 
the dark matter halo asymptotic circular velocity, V c , h, core 
radius r c ,h, an d 3 dimensional flattening, qs.h- The me- 
dian, 16th and 84th percentiles of the marginalized posterior 
PDFs for these parameters individually are given in Table [5] 

The constraints on stellar mass, dark halo density, and 
dark halo shape are discussed in more detail below. We first 
point out some of the main features of these figures. 



With strong lensing alone (Figure [TO]) , the halo param- 
eters are poorly constrained: The PDF for the core radius 
is almost uniform, while the PDFs for the halo velocity and 
halo axis ratio follow the priors. This is expected owing to 
the limited range in projected radius probed by the lens- 
ing constraints. There is, however, a good constraint on the 
stellar mass: log 10 (M»/M Q ) = 11.05t£;2§. This is a result of 
the axis ratio of the projected mass being quite low ( H5.2p . 

With kinematics alone (Figure [TTj) . the halo core radius 
is slightly better constrained, and the stellar mass is less well 
constrained. There is a strong degeneracy between the halo 
core radius and the stellar mass, with higher stellar masses 
requiring higher core radii — this is the classic disk-halo 
degeneracy. Related to this there is a degeneracy between 
halo velocity and core radius, with higher halo velocities 
requiring larger core radii. 
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Figure 11. Marginalised two-dimensional posterior PDFs for unconstrained mass model parameters (see Table[3]l using constraints from 
kinematics alone. In the histogram panels the vertical red lines show the median value. The median value together with the offsets to 
the 84th and 16th percentiles of the distribution is given in the top right corner. The priors are shown withe solid red lines. 



Adding the strong lensing constraints to the kinematics 
constraints breaks some of the degeneracies. Specifically, it 
removes the highest stellar mass solutions from the kinemat- 
ics only analysis. All posteriors are considerably tighter than 
the priors, illustrating the power of the combined analysis: 
for example, circular velocity is now known to 6% precision, 
and core radius is well constrained to be smaller than 5 kpc. 
However, there is still a degeneracy between halo core radius 
and stellar mass. There is also a residual degeneracy between 
stellar mass and halo shape — with low stellar mass solu- 
tions favoring oblate dark matter haloes. This degeneracy is 
expected as the total mass needs to be flattened to repro- 
duce the strong lensing (Section 15. 2\ . The flattening can be 
achieved with either a significant stellar disk component and 
a spherical halo, or a less massive disk and a more flattened 
halo. 



In Figure [T3] we show the rotation curves and strong 
lensing image predicted by two example mass models drawn 
from the posterior PDF given both lensing and kinematics 
data. These models both predict four images of the lensed 
source, including a faint counter-image that is consistent 
with the noise in the centre of the subtracted image. Both 
models' predicted rotation curves fit the outer part of the 
observed rotation curve very well; they also match well the 
central part of the observed rotation curve, which was not 
used in the fit. The two models have either the posterior 
median stellar mass, or much lower stellar mass; they can 
only be distinguished in the plot of the intrinsic, pre beam- 
smeared rotation curves, where the high stellar mass model 
has a significantly higher rotation velocity at radii less than 
one arcsec. This region could be probed with higher spatial 
resolution spectroscopy, or by making use of the velocity 
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Figure 12. Marginalised two-dimensional posterior PDFs for unconstrained mass model parameters (see Table[3]l using constraints from 
both kinematics and strong lensing. In the histogram panels the vertical red lines show the median value. The median value together 
with the offsets to the 84th and 16th percentiles of the distribution is given in the top right corner. The priors are shown with solid red 
lines. 



Table 5. Summary of fitted parameters: stellar mass (M 4 ); halo asymptotic circular 
velocity (v c h ); halo core radius (r c h); and 3D halo axis ratio (q 3 h ). 



log lo (M*/M ) Vch/Ikms" 1 ] r c , h /[kpc] q 3)h 
Lensing 11.05±g:§| 278l|J ^tgl °- 97 -o.2°s 



Kinematics 11.20t°'°l 275 ±33 ^-t'/r 1 -® 2 -t'%t 

Lensing + Kinematics 10.99t°;2g 276 ±i8 2 - 35 -i'.5 °- 91 iai3 
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Figure 13. Example mass model that fits the lensing (left panel) and kinematics (right panel) data. In the right panel two models are 
shown: median stellar mass (solid lines); low stellar mass (dashed lines). The red lines show the intrinsic model circular velocity, while 
the black lines show the model circular velocity after beam-smearing, finite slit width, and inclination effects are taken into account. 
Only the black points beyond 2 arcsec are included in the fit. The red point at small radii is the constraint from the stellar velocity 
dispersion, and disfavors the low stellar mass solution. 



dispersion information. Indeed, our cross-check point from 
Section [6. 51 would favour the high stellar mass model. 

In Figure [14] we show the inferred circular velocity pro- 
file, decomposed into baryonic and dark matter components. 
These estimates are based on the posterior samples using 
the joint lensing plus kinematics analysis. The solid lines 
show the median model from the posterior PDF, while the 
shaded regions enclose 68% of the models. In the radial 
region where we have observational constraints (i.e., from 
the Einstein radius to the last rotation curve point) the 
total circular velocity is well constrained. For example, at 
2.2 disk scale lengths (8.1 kpc), the total circular velocity 
is V2.2 = 289 ± 4kms _1 . The circular velocity profiles of 
the baryons and dark matter are not as tightly constrained. 
Nevertheless, we can still infer interesting constraints on 
the dark matter fraction as a function of radius, and thus 
determine whether or not SDSS J2141— 0001 has a maxi- 
mum disk. A working definition of a maximum disk is when 
V r di S k(2.2 J Rd)/Kot(2.2 J R d ) = 0.85 ±0.10 (Sackett 1997). Here 
Vdisk(2.2i?d) is the circular velocity of the disk at 2.2 disk 
scale lengths, and Kot(2.2i?d) = V2.2 is the total circular 
velocity at 2.2 disk scale lengths. 

A galaxy may have a sub-maximal disk, but still have a 
maximal baryonic component due to the bulge. Thus we con- 
sider the contribution of the baryons (i.e., bulge plus disk) 
to V2.2 to be of more relevance than just the contribution of 
the disk to V2.2. We find that Vb a r(2.2 J R d )/V r 2 .2 = 0.67t l°, 
which suggests that SDSS J2141— 0001 is sub- maximal at 2.2 
disk scale lengths. However, the baryon contribution to the 
total circular velocity increases towards smaller radii (lower 



panel of Fig. \U\) such that V ba *{Rb) /Vtot = 0.99l°;o9, and 
thus SDSS J2141-0001 is maximal at the bulge half-light 
radius. Converting circular velocities into spherical masses, 
results in a dark matter fraction of /dm = 0.55t 'i5 within 
2.2 disk scale lengths, and /dm = 0.02^0 ^2 within the bulge 
half-light radius. 

7.2 Constraints on the stellar IMF 

Figure [TS] shows the posterior PDFs from our joint lens- 
ing and kinematics analysis together with those from SPS 
models for both Chabrier and Salpeter IMFs. From our lens- 
ing and kinematics analysis the stellar mass of the galaxy 
is found to be log 10 M„ = 10.99lo'25- This is in excellent 
agreement with the stellar mass derived from SED fitting 
assuming a Chabrier IMF, which is log 10 M* = 10.97T0.07. 
A Salpeter IMF results in stellar masses 0.24 dex higher. Our 
analysis thus mildly favors a Chabrier IMF over a Salpeter 
IMF. We can quantify this agreement by integrating the 
likelihood over a prior for the stellar mass defined by either 
the Chabrier or the Salpeter SPS model PDFs. The ratio of 
these integrals is the Bayes factor, or evidence, in favour of 
a Chabrier IMF; we find its value to be 2.7, which is to say 
that the data are 2.7 times more likely to have come from 
the Chabrier model than from the Salpeter one. If these are 
the only two models possible, then there is a 73% chance 
that the Chabrier model is the true one. This corresponds 
to weak evidence in favor of Chabrier vs Salpeter. 

In our current analysis we ignore the possibility of cold 
gas. For massive spiral galaxies the cold gas fraction is ~ 



© 2011 RAS, MNRAS 000,[TJf23] 



18 Button et al. 



300 HS b 

.1 : % 



200 



>° 100 L 



> 



> 




I ^-1 I I i I I I I L I I I I I I 

5 10 15 

Rodius [kpc] 

Figure 14. Circular velocity profiles (upper panel) and spherical baryon fractions (lower panel) from our joint lensing plus kinematics 
analysis. The solid lines show the median, while the shaded regions enclose 68% of the posterior PDF. The total circular velocity (black 
line and grey shaded region) is well constrained outside of the Einstein radius, iiEin, and up to the last rotation curve point at 13 kpc. 
The contributions of the baryons (red lines and shaded regions) and the dark matter (blue lines and shaded regions) are more uncertain. 
However, at the bulge half-light radius, Ri,, the galaxy is baryon dominated (and thus is "maximal"), while at 2.2 disk scale lengths the 
baryons fraction is roughly 50% (and thus is "sub-maximal"). 
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Figure 15. Inference on stellar mass from lensing and kinematics 
(histogram) compared with SPS models (solid lines) assuming a 
Chabrier IMF (black) and Salpeter IMF (red). The Bayes factor 
in favor of a Chabrier IMF, compared to a Salpeter IMF is 2.5. 



20 ± 10% (assuming a Chabrier IMF), split roughly equally 
between atomic and molecular gas (e.g., Dutton & van den 
Bosch 2009). If the cold gas is distributed like the stars, then 
the lensing+kinematics stellar mass is actually a baryonic 
mass, greater than or equal to the actual stellar mass. If 
the cold gas is more extended than the stars, as is often the 
case, then we will still be over-estimating the stellar mass, 
but by a smaller amount. To estimate an upper limit to the 
impact of cold gas on our derived stellar masses we assume 
that the gas mass for SDSS J2141— 0001 is distributed like 
the stars. For each model in the posterior PDF we draw a 
gas mass from a log-normal distribution centered on M gaa = 
1.8 x 1O 1O M , with a standard deviation of 0.3 dex. We then 
subtract off the gas mass from the gas free stellar mass to 
derive the "true" stellar mass. The results of this exercise are 
shown in Fig. 1161 The resulting median and 68% confidence 
interval on the stellar mass is log lo (M*/M0) = 10. 891° 331 
i.e., 0.1 dex lower than when ignoring the cold gas. The 
Bayes factor in favor of a Chabrier IMF over a Salpeter 
IMF has increased from 2.7 to 11.9, which corresponds to 
strong evidence. Thus by ignoring the cold gas we could be 
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Figure 16. Effect of gas mass on the inferred stellar mass from 
lensing and kinematics. For each model galaxy in the posterior 
PDF we draw a gas mass from a log-normal distribution with 
mean and standard deviation typical for massive spiral galaxies. 
The resulting PDFs for the stellar and gas mass are shown as 
black hatched histograms. For comparison, the blue shaded his- 
togram shows the posterior PDF on the stellar mass assuming no 
gas mass. Thus accounting for cold gas mass reduces the stellar 
mass derived from lensing and kinematics by ~ 0.1 dex. 



over estimating the stellar mass by ~ 0.1 ± 0.05 dex, which 
strengthens the case for an IMF lighter than Salpeter. 

How does this result compare to previous work? Us- 
ing maximal disk fits to spiral galaxy rotation curves in the 
Ursa Major cluster, Bell & de Jong (2001) placed an upper 
limit on the stellar mass-to-light ratio normalisation, favor- 
ing IMFs with stellar masses 0.15 dex lower than Salpeter, 
the so-called diet-Salpeter IMF. We note that a Salpeter 
IMF is also disfavored for fast rotating elliptical galaxies 
(Cappellari et al. 2006; Treu et al. 2010; Auger et al. 2010; 
Barnabe et al. 2010), but is favored for massive elliptical 
galaxies (Treu et al. 2010; Auger et al. 2010; van Dokkum 
& Conroy 2010). Thus comparing our result with those for 
massive ellipticals, supports the idea that the IMF is not uni- 
versal, but dependent on galaxy mass and/or Hubble type. 

By shifting our Salpeter stellar mass PDF by —0.15 dex, 
we find that for SDSS J2141-0001 a diet Salpeter IMF cor- 
responds to a stellar mass of log 10 (M*/MQ) = 11.06 ± 0.07. 
This IMF is favored over a Salpeter IMF, by Bayes factors 
of 3.5 (assuming no cold gas), and 9.9 (assuming a gas mass 
of log lo (Af gas /M ) = 10.26 ± 0.30). There is little evidence 
distinguishing between diet Salpeter and Chabrier IMFs. 



7.3 Constraints on the dark halo density profile 

N-body simulations have shown that in ACDM cosmologies 
dark matter haloes are expected (in the absence of baryonic 
effects) to have very specific structure. The mass density 
profiles can be well approximated by the so called NFW 
profile (Navarro et al. 1997). This has a density profile that 



varies from p(r) oc r _1 at small radii, to p oc r~ 3 at large 
radii. The radius where the logarithmic slope of the density 
profile is dlnp/dlnr = —2 is known as the scale radius, r s . 
The halo scale radii are tightly correlated with the virial 
masses of dark matter haloes, M V1I . This correlation is usu- 
ally expressed in terms of the halo concentration, c = r v i T /r s , 
where r v i r is the virial radius. Halo concentrations are only 
weakly dependent on halo mass, with a relation of the form 
c oc M~® A (Maccio et al. 2007). The scatter in halo concen- 
tration, at fixed halo mass, for relaxed haloes is small ~ 0.11 
dex (Jing 2000; Wechsler et al. 2002; Maccio et al. 2007). 

Observationally, measuring halo concentrations is a 
challenge because halo virial masses are poorly constrained 
for individual galaxies, due to the lack or sparsity of dynam- 
ical tracers at large radii. A more observationally accessible 
measure of dark halo structure is the parameter Ay/ 2 , which 
depends on the maximum halo circular velocity, Umax, and 
the radius where circular velocity of the halo is half of the 
maximum, r v / 2 (Alam, Bullock & Weinberg 2002): 

(22) 
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= 5 x 10 a 



For NFW haloes there is 
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ry/ 2 / [/i _1 kpc] 

a one-to-one mapping between 
and thus one can compare the 
observed A V / 2 with predictions for ACDM haloes. Figure lTTl 
shows the predictions for A v / 2 — U max in a WMAP 5th year 
cosmology (Dunkley et al. 2009) from Maccio, Dutton, & van 
den Bosch (2008) . The shaded regions show the 1 and 2a in- 
trinsic scatter. The large symbols show measurements from 
dwarf and low surface brightness galaxies after subtracting 
of the baryons (de Blok et al. 2001; de Blok & Bosma 2002; 
Swaters et al. 2003). These are in excellent agreement with 
the predictions from ACDM. 

In our mass model of SDSS J2141— 0001, the dark mat- 
ter halo has a softened isothermal density profile, which has 
U m a X = U c ,h, and r v / 2 = 1.1263 r Cjh . We can compute A v / 2 
for this model and compare it with the NFW profile ha- 
los of the simulations. For our model the median and un- 
certainty (corresponding to 16th and 84th percentiles) is 
log 10 A v / 2 — 5.91q 5- The median is 2.7cr higher (in terms 
of intrinsic scatter) than that predicted for pristine ACDM 
haloes, although the full posterior PDF overlaps the ACDM 
predictions, as shown in Figure 1171 The 16th percentile of 
the PDF for A v / 2 only corresponds to a Tier deviation from 
the ACDM distribution. Thus there is a suggestion that the 
SDSS J2141-0001 halo is higher density than expected. We 
note that, as shown in Fig. 1181 the central density is highly 
correlated with the stellar mass, with lower A v / 2 for higher 
M*. 

There are two interpretations of the higher than ex- 
pected halo density. 1) The halo has undergone contrac- 
tion in response to galaxy formation (e.g., Blumenthal 
et al. 1986). 2) The halo of SDSS J2141-0001 has a higher 
density than typical haloes of the same mass. In order to 
distinguish between these two scenarios it is necessary to un- 
derstand the selection function of the SWELLS lens galax- 
ies. We cannot do this with one galaxy, and therefore we 
will leave this for future work. Nevertheless, we can gain 
some insight by investigating where SDSS J2141— 0001 falls 
on disk galaxy scaling relations. 

We consider the relations between stellar mass (Af»), ro- 
tation velocity at 2.2 disk scale lengths (U2.2), and disk scale 
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Figure 17. Central density of the dark matter halo, Ay/ 2 vs 
maximum halo circular velocity, Vmax- The solid black line shows 
the prediction for pristine dark matter haloes in the concordance 
ACDM cosmology (WMAP5) from Maccio ct al. (2008). The 
shaded regions show the 1 and 2<r intrinsic scatter. The colored 
symbols show measurements from dwarf and low surface bright- 
ness galaxies, after subtraction of the baryons (de Blok et al. 2001, 
dBOl; de Blok & Bosma 2002, dB02; Swaters et al. 2003, S03). 
The small black dots show samples from the posterior PDF for 
SDSS J2141-0001. The large red dot and error bars show the 
median and 68% ranges of this PDF, respectively. 

length (R d ). For SDSS J2141-0001 the values of these pa- 
rameters are V2.2 — 289 km s -1 , log 10 (M*,ciiab/M©) ~ 11.0, 
and Rd — 3.7kpc. For a stellar mass of lO^M© we ex- 
pect V2.2 = 229±25kms~ 1 (Dutton et al. 2010b). Alterna- 
tively, for a rotation velocity of V2.2 = 289 kins -1 we expect 
a stellar mass of log 10 (M t /M Q ) = 11.39lo'ig, assuming a 
Chabrier IMF. For a rotation velocity of V2.2 = 289 km s" 1 
we expect Ra = 6.0^'gkpc (Courteau et al. 2007; Dutton 
et al. 2007). Thus SDSS J2141-0001 is offset to low stellar 
mass (by ~ 2a and small size (by ~ 1.3a) at fixed rotation 
velocity, which means that it has a higher baryonic and total 
mass density than typical massive spiral galaxies. 

The high central density of SDSS J2141-0001 may be 
related to how it was selected. If the central densities of spi- 
ral galaxies are close to the critical value for strong lensing, 
then the galaxies that are observed to be strong lenses could 
be a biased sub-set of the population. Thus the conclusions 
we draw for SDSS J2141— 0001 might not necessarily be ap- 
plicable to spiral galaxies in general. This selection effect is 
likely to affect the interpretation of some parameters more 
than others. For example, in order to interpret the central 
densities of dark matter haloes in terms of the halo response 
to galaxy formation it is necessary for the lenses to reside in 
an unbiased subset of haloes. However, if we assume that the 
IMF is universal, at least for massive spiral galaxies, then the 
constraints we place on the IMF using spiral galaxy strong 
lenses will be independent of any selection bias. Larger sam- 
ples of spiral galaxy lenses are obviously needed to fully 
characterise the selection biases in spiral galaxy lenses. 



Figure 18. Central density of the dark matter halo, A^/ 2 vs 
stellar mass, M„. The solid black line shows the prediction for 
pristine dark matter haloes in the concordance ACDM cosmol- 
ogy (WMAP5). The shaded regions show the 1 and 2<r intrinsic 
scatter. The black dots show samples from the posterior PDF for 
SDSS J2141— 0001. The maximum likelihood stellar masses from 
SPS models assuming Chabrier, diet-Salpeter, and Salpetcr IMFs 
are shown with vertical dashed lines. 



7.4 Constraints on the dark halo shape 

N-body simulations have shown that in ACDM cosmologies 
dark matter haloes are triaxial, with a preference towards 
prolate shapes (Jing & Suto 2002; Bailin & Steinmetz 2005; 
Kasun & Evrard 2005; Allgood et al. 2006; Bett et al. 2007; 
Maccio et al. 2007, 2008). For a halo mass of M vir = 1O 12 M , 
which is typical for massive spiral galaxies, the minor to ma- 
jor axis ratio ratio c/a ~ 0.63 ±0.1, and the intermediate to 
major axis ratio b/a ~ 0.80T0T (Maccio et al. 2008). ACDM 
haloes are also found to be more prolate at smaller radii 
(Allgood et al. 2006; Maccio et al. 2008; Abadi et al. 2010) . 

The assembly of a central galaxy is expected to modify 
the three dimensional shape of the dark matter halo (Katz & 
Gunn 1991; Dubinski 1994). Using cosmological simulations 
Abadi et al. (2010) found that as a result of the assembly 
of the central galaxy the haloes become nearly oblate with 
b/a ~ 0.95 and c/a ~ 0.85, and the axial ratios become ap- 
proximately independent of radius. Similar results have been 
obtained from other cosmological simulations (e.g., Tissera 
et al. 2010). It should be noted that these simulations over- 
predict the baryon to dark halo mass ratios, and thus the 
effect of galaxy assembly on the halo shapes could be over- 
estimated. 

In our mass models we assume the halo is axisymmet- 
ric, with 3D axis ratio qa,h- A spherical halo has 53^ = 1, 
an oblate halo has q^^ < 1, and a prolate halo has 53^ > 1. 
We adopt a log-normal prior on the halo axis ratio, cen- 
tered on 53. h = 1. For our fits to lensing only or kinematics 
only the posterior PDF for q^^ is identical to the prior, but 
for the joint analysis slightly oblate haloes are preferred: 
g3,h = 0.9llo'}3, and prolate haloes with q3,h > 1.2 are 
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strongly disfavored. Thus our results for SDSS J2141-0001 
support the notion that galaxy assembly sphericalizes the 
dark matter halo, and perhaps even flattens it towards the 
disk. 



8 CONCLUSIONS 

We have presented an analysis of the strong gravitational 
lens SDSS J2141-0001, discovered as part of the SLACS 
survey, using data from HST and the Keck telescopes. The 
lens galaxy is a high inclination disk dominated galaxy with 
Tf'-band bulge fraction of 0.2, showing stellar rotation in 
multiple spectral lines. A singular isothermal ellipsoid lens 
model provides a circular velocity of V c = 254+]^ km s _1 
and an axis ratio of q = 0.421q[ 12 . 

We perform a joint fit to the multi-filter surface bright- 
ness, lensing and kinematics data using a self-consistent 3- 
component mass model, and from it draw the following con- 
clusions: 

• The lensing and kinematics constraints yield a stellar 
mass of log 1Q (M,/M ) = 10.99±g;|^ (68% confidence inter- 
val), independent of the IMF. 

• This value is in excellent agreement with the stellar 
mass derived from the SED using SPS models and assum- 
ing a Chabrier (2003) IMF: Iog 10 (M*/M Q ) = 10.97tg;J£. 
A Salpeter (1955) IMF results in stellar masses 0.24 dex 
higher: our analysis marginally favors a Chabrier IMF over 
a Salpeter IMF, by a Bayes factor of 2.7. 

• Accounting for the expected gas mass reduces the lens- 
ing and kinematics stellar mass by 0.10 ± 0.05 dex, and in- 
creases the Bayes factor in favor of a Chabrier IMF to 11.9. 

• At 2.2 disk scale lengths the spherical dark matter frac- 
tion is /dm = 0.55^0 i5i suggesting that the baryons are 
sub- maximal. 

• The dark matter halo has a maximum circular veloc- 
ity of V c ,h = 276tig kms -1 , an d a core radius of r Cj h = 
2.4l^'jkpc. The corresponding central density parameter 
log 10 A v /2 = 5.9^0 5 is higher than expected for uncon- 
tracted NFW haloes in the concordance ACDM cosmology, 
which have log 10 A v / 2 = 5.2 and an intrinsic scatter of 0.3. 

• This high density could either be evidence for halo con- 
traction in response to galaxy formation (e.g., Blumenthal 
et al. 1986), or the result of a selection bias towards high con- 
centration haloes. A larger sample with well-characterised 
selection function is required to make further progress. 

• The dark matter halo is oblate, (/3,h = 0.91 jlg'^, with 
a probability of 69%. This finding provides support for the 
notion that galaxy assembly turns strongly prolate triaxial 
dark matter haloes into roughly oblate axisymmetric haloes 
(e.g., Abadi et al. 2010). 
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APPENDIX A: THE CHAMELEON 
APPROXIMATION TO A SERSIC PROFILE 

In this appendix we derive an approximation to a Sersic pro- 
file as the difference of two non-singular isothermal ellipsoids 
(NIE's). 

The Sersic profile is specified by three parameters: a 
normalization, a radial scale, and a shape parameter com- 
monly known as the Sersic index, n. In its simplest form it 
is given by: 



E(7?) = E exp 



l/r. 



(Al) 



where Eo is the central surface density, and Ro is the radial 
scale. The Sersic profile is commonly expressed in terms of 
the effective radius, R B , which encloses half of the projected 
mass, and the effective density, E c = E(J? C ): 



E(i?) = E c exp <^ -h 



(A2) 



Where the term b„ » 2n — 0.32. Here we use the asymptotic 
expansion of Ciotti & Bertin (f999) to 0(n~ a ) valid to one 
part in ~ 10 4 for n > 0.36: 



b„ = 2n 



1 



405™ 



+ 



46 



+ 



131 



25515n 2 1148175n 3 
2194697 . _ 5 , 



+ 0(n" 



(A3) 



30690717750n 4 

Equation I All and Equation I A2I are related to each other via 
R c = (b n ) n Ro (A4) 

and 



E c = exp(— 6„)E . 



(A5) 



The Chameleon profile is the difference between two 
NIEs with different core radii, but the same normalization 
(which insures the total mass is finite): 



Echm (R', El , R\ , R2) — 

Ei Ei 



^R? + R\ ^R 2 +R 



(A6) 



For the purpose of comparing to a Sersic profile, we use the 
following parametrization. 



Eo 



Ro 



E c hm (J?; Eo, J?o, a) = 
Ro 



1 



y/R? 



R 2 ^/R 2 + (Ro/a) 



(A7) 



Where Eo = Ei(l — a)/R\ is the central surface density of 
the chameleon profile, Ro = Ri, and a is the ratio between 
the core radii of the two NIE's: a — R1/R2, so that < a < 
1. The total mass of the Chameleon profile is 



Mciim = 7rE -R / a - 



(A8) 



We wish to find the correspondence between the Sersic 
parameters, Eo,sor, Ro.Sor, ^Sor, and the Chameleon pa- 
rameters, E ,chm, Ro.Chm, QChm- We do this by fitting a 
Chameleon profile to a Sersic profile. We are interested in 
mass profiles, so we fit R E(i?) to give appropriate weight to 
the density profile at large radii. Figure [ATI shows the best 
fit parameters as a function of Sersic index. We fit these 




Figure Al. Relation between Chameleon and Sersic parameters 
as a function of Sersic index. 



Table Al. Parameters of cubic fitting formula (Equation I A9|| to 
the relations in Figure [ATI 



,'/() 



Vi 



U2 



lo sio 
l°gio 



S 0,Chn 
E 0,Scr 



-Ro.Sor 
lo gl0("Chm) 



1.69 -0.254 

1.15 0.078 
2.03 -0.739 



-0.259 
-0.184 
-0.527 



-0.036 0.014 
0.473 -0.079 
-0.012 -0.008 



relations between chameleon and Sersic parameters with a 
cubic function: 



y = yo 



+ yi(x- xo) + 2/2(2: - x ) 2 + 1/3(2; - x ) 3 (A9) 

The parameters of these fits are given in Table IA1I Our 
fitting function is valid for 1 ^ ns C r ^ 4. 

Figure [A"2l shows Chameleon fits to Sersic profiles with 
Sersic index n = 1,2,3, & 4 using the fitting function in 
Table IAT1 The left panels show log 10 E(i?) vs (R/R c ) 1/n . In 
these units Sersic profiles are straight lines. The Chameleon 
profiles deviate most significantly from Sersic profiles at 
small radii, where the Chameleon profile has a constant den- 
sity core. The middle panels show RS(R) vs R/R e , which 
are the curves that were fitted against. The right panels show 
the cumulative mass profile. For radii between 0.5 < R < 3.0 
effective radii the mass residuals are only a few percent, and 
thus for most strong lensing and dynamical purposes the 
Chameleon approximation to a Sersic profile is of sufficient 
accuracy. 

This paper has been typeset from a TpX/ F/IpX file prepared 
by the author. 
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Figure A2. Chameleon fits (red solid lines) to Sersic profiles (dotted lines) with n = 1, 2, 3, 4. The Chameleon profile reproduces the 
cumulative mass profile of a Sersic profile to a few percent over most radii interesting for strong lensing and kinematics. 
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